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The scale dependent intermittency exponents in developed 
hydrodynamic turbulence are calculated assuming a natural 
hierarchy of correlations in the turbulence. The major corre- 
lations are taken into account explicitly, while the remaining 
small correlations are considered as perturbations. The re- 
sults agree very well with the currently available experimental 
data. 

PACS-number: 47.27 

• INTRODUCTION: The energy dissipation in 
developed hydrodynamic turbulence is strongly non- 
uniform in space (and time). This property is usually 
referred to as an intermittency. A reason for intermit- 
tency was already present in the original Kolmogorov- 
Obukhov model of turbulence, as it was uncovered by 
the well-known comment of Landau pi. The major pro- 
cess in developed turbulence is the energy transfer from 
turbulent motions of large spatial scales (where the en- 
ergy is basically contained) to the smallest scales (where 
the dissipation occurs). The transferred energy flux be- 
comes more and more spatially non-uniform as it goes 
to smaller and smaller scales, which ultimately leads to 
pronounced intermittency. A reason for this is that the 
velocity field at smaller scales is capable of changing in 
time much faster than at larger scales. It implies that 
a quasi-stationary (in statistical sense) picture of turbu- 
lence has time enough to be established at smaller scales 
before a noticeable change in the local large-scale flow oc- 
curs. Therefore, the statistical characteristics of smaller- 
scale motions depend just on the local (in space and time) 
larger-scale motion. Then, there is no smoothing mech- 
anism for spatial fluctuations in the density of energy 
flux and those are amplified as the flux goes from larger 
to smaller scales. The quantitative description for such 
an intermittent amplification has not been developed yet 
because of the complexity of the problem. 

An idea is suggested below on how to overcome the dif- 
ficulties, assuming that not all of the correlations causing 
the intermittency are equally important. It appears to 
be possible to derive and solve explicitly an equation that 
takes into account the major correlations. After a renor- 
malization, a local flux-conservation constrain leads to an 
equation that has (in addition the Kolmogorov solution) 
an intermittent solution. The latter agree very well with 
experiment, which supports the employed assumptions 
about hierarchy of correlations. 

• ENERGY FLUX AND DISSIPATION: The 



problem of intermittency can be considered quantita- 
tively in the terms of energy and flux densities in (fc, r)- 
space. Provided the energy is pumped at the smallest 
wave numbers fc ~ feo, damped at the largest k ~ kd ^$> fco 
and there is neither pumping nor damping at all inter- 
mediate fc, the densities of energy and flux satisfy in the 
latter - inertial - range the continuity equation. It can 
be averaged over directions of k and over some domain of 
fixed shape in r-space. The result contains a surface in- 
tegral (over the boundary of domain) , that describes the 
spatial transfer of energy from the domain to its neigh- 
borhood, and an integral over the domain, that repre- 
sents the averaged density of energy flux from smaller 
to larger k. The latter quantity is denoted below as 
e(R, x, fc), where R and x define the spatial position and 
linear size of the domain; the time variable and variables 
describing the shape of domain are omitted for simplicity. 
One could expect that the average flux density e(i?, x, k) 
tends at k 3> 1/x (but k is inside the inertial range) to a 
fc- independent limit. Such a limit, if it exists, should be 
identified with the average density of energy dissipation 
in given domain at given moment of time, which quantity 
is denoted below as e(R,x). 

Statistical properties of random fields e(R, x, fc) and 
e(R,x) may be sensitive to the shape of domain used 
for the averaging of energy flux and dissipation densities. 
Therefore, for comparison with experiments, the same 
shapes as there must be employed in the theory as well. 
Because of the technical reasons, the experiments usually 
deal with the so called "one-dimensional cuts" , obtained 
by means of small probes that record the longitudinal 
velocity of fluid in fixed points passed by the stream. In 
some recent experiments the transverse component of ve- 
locity is measured, but the currently available database 
contains primarily records of the longitudinal velocity. 
When such data are processed, the density of energy dis- 
sipation rate is substituted usually by the square of the 
velocity derivative along the record. 

• BASIC EQUATIONS: In order to describe the 
amplification of fluctuations in density of energy flux 
from larger to smaller scales, consider the following con- 
ditional probabilities. 

Let G n (M)dM be the probability to find the flux 
e(R,x/2 n ,2 n /x) in the range (M,M + dM)e, provided 
e(R, x, 1/x) = e. The ensemble averaging means here and 
further the i?-space averaging. The probability density 



G n (M) does not depend on x (inside the inertial range). 
It is a scale- invariance hypothesis. The e-dependence of 
G n (M) (which would imply a correlation between the 
local amplification coefficient M and energy flux den- 
sity e) is assumed to be negligible in a zero-order ap- 
proximation. When correlation between the consequent 
amplification coefficients is also negligible, the moments 
G n . q = J dMM q G n {M) depend on n exponentially: 
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Let g n (M)dM be the probability to find the flux 
e(R,x,2 n /x) in the range (M,M + dM)e, provided 
e(R,x,l/x) = e. The consequent amplification coeffi- 
cients of this kind are strongly correlated, so that n- 
dependence of g n {M) is not trivial. There is, however, 
a relation between the probability densities g n (M) and 
Gi(M) = G(M). To find it, consider the "one dimen- 
sional cut" of length x as the sum of its halves (which 
length is x/2). The conditional probability to find the 
flux e(R, x/2, 2/x) in the range (M, M + dM)e, provided 
e(R,x, l/x) = e, is G{M)dM. It relates to each half of 
the segment x, when the dependence on the position of 
the half inside the segment is negligible. The amplifica- 
tion coefficients for the halves virtually do not affect each 
other, because of the strongly fluctuating short-scale spa- 
tial transfer of the energy in directions transverse to the 
segment. Then, the following recurrent relation holds up 
to small corrections: 

g n (M) = / dMidM 2 dM 3 dM 4 G(Mi)G(M 2 ) x 
g n -i(M a )g n _ 1 (M A )8 (M - {M X M Z + M 2 M 4 )/2) . (2) 

It can be reduced to the recurrent relation for moments 
g, hq = / dMM q g n (M), that looks as 



<r— s q\G qi G q - qi g n -i, qi g n -i 
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Eq.(g) with "initial conditions" go. q = g n fl = 1 allows 
one to express g n ^ q in the terms of G q for all positive 
integer n and q. 

The n-dependent moments of energy flux densities av- 
eraged over segment x and over its part x/2 n at the same 
wave number k n = 2™ jx would be equal to each other for 
a spatially uniform flux. For an intermittent energy flux, 
the ratio 



^n,q / 9w 



2 n ^n,q 
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increases in n (which indicates the amplification of fluc- 
tuations). According to (||) and (||), the "intermittency 
exponents" [i. n ,q are given by 



Mn,g 



log 2 G q ~ n l log 2 g n 



(5) 



• EXPLICIT SOLUTION OF EQ.@: Equation 
S) can be solved explicitly, and thus the "intermittency 



exponents" /j, n>g can be expressed in the terms of n- 
independent functions G q . Provided the latter satisfy 
restriction G q < 2 9 ~ 1 , the solution g n ^ q tends to finite 
limit g q as n tends to infinity. Then, the "intermittency 
exponents" /i„ j( j tend to /x g = log 2 G q as n tends to in- 
finity, and the above restriction on G q is equivalent to 
H q < q — 1. Explicit formulas for g n>q that solve (0) can 
be derived consequently for q = 1, 2, .... 

For q = 1, taking into account that Gq = G\ = 1, one 
gets from (0) g n ^i — 1, and hence //„i = 0. 

For q = 2 the solution of (0) is 



9n.2 = g-2 + c 2 (G 2 /2) ; 
g 2 = 1 1 (2 - G 2 ) , c 2 = 1 - g 2 . 



(6) 



For q = 3 the solution is 

9n,3 =93 + c 3 (G 3 /4f + c 32 (G 2 /2) n ; 
9a = T^f, c 32 = ^ G ~ 3 , c 3 = l- 33 - C32 . (7) 

For q = 4 it is 

9n,A = 94 + C 4 (G4/8)" + C 43 (G 3 /4)" + 

+ c 42 (G 2 /2) n + c 422 (G 2 /2) 2 "; 

4G393+3G|s| 4c 3 G 3 

94 — 8-G4 ' 43 2G3-G4 ' 

_ 4c 3 2G 3 +6c 2 Glg2 3c^ 

C 42 - 4G 2 -G 4 ' C 422 - 2G 'i Gi , 



'2- Lj4 

c 4 = f - <? 4 - c 43 - c 42 - c 422 



(8) 



For g 



9n,5 = ff5 + c 5 (G5/I6)" + c 54 (G 4 /8) n + c 53 (G 3 /4)' 
+ c 52 (G 2 /2)" + c 532 (G 2 G 3 /8)" + c 522 (G 2 /2) 2n ; 

„ _ g G4g4+2G 2 G 3 g2g3 
55 - ° ie^Gl ' 



■■ !■ ■ i ,,„„ r; C43G.i+2c3G3G2g2 

C54 _ 2Gi-G 5 ' C53 ~~ ° 4G3-G5 
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C42G 4 +2c2G3G2g3+2c32G 3 G2g2 



^4 — l_r5 

C52 = 5 



8G2 — G5 



C 532 - 2 G 2 G3-G 5 ' C522 - b 4G*-G 5 



IOC2C3G2G3 
2G2G3 — G5 

C5 = 1 - 55 - c 54 - c 53 - c 52 - c 532 - c 522 . 



(9) 



The similar (but longer) formulas exist for q = 
6,7,8,.... In particular, the n-independent function g q 
is given in general by 



9i 
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Consider now the alternative option G q > 2 q ~ 1 . Un- 
der such a condition, g nq does not tend to a finite limit 
as n tends to infinity, but increases exponentially like 
9n, q oc {G q /2 q - 1 ) n . It follows then from (|) that the in- 
termittency exponent \x n ,q tends to q — 1 as n tends to 
infinity (and one can see that ix n ,q approaches to q — 1 
from below). There is a simple physical interpretation 



of such a behavior. It indicates that the strongest fluc- 
tuations of e{R,x) increase faster than 1/x as x tends 
to zero (so that the density of energy dissipation rate 
has singularities non-integrable over the one-dimensional 

cuts). According to the estimate e(R,x) ~ v(R,x) jx 
(where v(R, x) is a smoothed velocity variation in a fluc- 
tuation of size x around point R), it implies that the 
sweep v(x) ofv(R, x) variation in i?-space increases when 
the fluctuation size x tends to zero. Thus, the alterna- 
tive G q > 2 q ~ 1 corresponds to blow-up of velocity itself, 
rather than just its gradients, in Euler fluid. 

In spite of many efforts applied to the problem, there is 
no clear theoretical answer yet to the question, whether 
velocity blows-up or not in Euler fluid. Practical an- 
swer can be extracted from the experimental data on 
intermittency, provided the collected statistics is suffi- 
cient to detect the strongest blow-up. The currently 
available database supports assumption [i q < q — 1 for 
q = 2,3,4,.... 

• DEGENERATE SOLUTIONS: The above for- 
mulas for g ntq contain all possible items of the kind 
(Gq/2?- 1 )™ , q = 2,3,... and their products. There is 
however a special class of solutions that contain only the 
items with q = 2. It corresponds to the physical situa- 
tion when only one major kind of correlations exists, that 
dominates over all other correlations. For such degener- 
ate solutions, all quantities G q (or fi q — log 2 G q ) can be 
expressed in the terms of fj,2- In the limit of infinitely 
small [ii-, the expression looks as 



M<? ~ 1*29(9- l)/2 



(11) 



The dimensionality D q = 3 + \i q — q\s! q corresponding to 
(O) is D q ps 3 — fj,2 Q 2 /2 . It would turn into zero at 
q = qu ~ (6/M2) 1 / 2 , if ( PI ) would be applicable there, - 
which is not so, as the applicability condition /i q <C 1 is 
violated at q ~ qu- In fact, the applicability condition 
is a little bit more soft quantitatively, it is G q — 1 <C 1. 
The slightly improved formula (|llj) is 



Gq-l 



(GV 



-1)9/2 



(12) 



It appears to be of a reasonable accuracy for real turbu- 
lence up to q = 5 (see more accurate formulas below for 
the comparison). 

• INTERMITTENCY EXPONENTS: For a 

moderately small y,i that occurs in real turbulence, the 
above formulas require a modification. An appropriate 
formula for fi q is derived from very general heuristic rea- 
sons in [0 (where the major mechanism of intermittency 
is supposed to be blow-up of velocity gradients in Euler 
fluid, leading to point singularities). This formula looks 
as 



M<? = 7 9 I arsn (<l/ a ) ~ arsn (V a )] > 1 < 1m 
Hg = a M q - 3 , (a M = Mg M )> 9 > 9m ; 

, MA 
9m = 3/7 



l/2+(l/4+(a 7 /3) 2 ) 



1/2 



(13) 
(14) 

(15) 



where arsh£ = ln(£ + -\/C 2 + 1) is the hyperbolic arcsin, 
prime signifies g-derivative, 7 and a are some numbers. 
For a~^> q, (13) reduces to ( |il| ) with ^ 2 = 27/a. 

There is a connection between parameters 7 and a, so 
that fiq can be expressed in the terms of \xi for realistic 
values of the latter parameter as well. The extra condi- 
tion has the form 



F(aM 7 ao) 



■)3-a 



7 = 0(1) 



(16) 



(which is easy to understand in the framework of a cubic 
shell model). The quantity a q = fi' can be referred to 
as the "singularity exponent" . When q increases from 
zero to <7mj ot q increases from its smallest value a^ < 
to its largest value aw = a qM . The exponent aM 
corresponds to the strongest possible amplification of e in 
a "hot spot", i.e., to the strongest blow-up in Euler fluid. 
The exponent ao < corresponds to the background 
that goes down to compensate for amplification of e in 
the hot spot and to secure the energy flux conservation. 
The flux conservation condition around the strongest hot 
spot gives 



F(atM,ato) 



A ■ 



(17) 



For comparison between the theory and experiment, 
the above definition of intermittency exponents /i„. g (se- 
lected to get explicit analytical formulas) must be slightly 
modified and also small effects neglected in the derivation 
of eq.(0) must be taken into account. The different defi- 
nitions of intermittency exponents can be identified with 
each other, as long as only scale-independent intermit- 
tency exponents \x q (corresponding to infinite scale ratio) 
are considered. The really measured intermittency expo- 
nents should be considered taking into account the finite 
scale-ratio effects which are sensitive to the definition. 
In particular, the above definition must be distinguished 
from that used in || . The latter deals with the amplifica- 
tion coefficients M(R, x, y) = e(R, y)/e(R, x), and defines 
intermittency exponents in a way equivalent to 



M(R,x,x/2 n )i = 2 n ^-" 



(18) 



Here bar signifies the R-sp&ce averaging. The result does 
not depend on x due to the scale invariance. 

The modification can be neglected for infinitely small 
/i2- For a moderately small fj,2, the quantities /i„ iQ defined 
by (n8|) are expressed in the terms of \xi by 



/<,,- — (log 2ff „, g ) C , Cwl+Ml-tfO (I 9 ) 



(where fi q is given by (|lq)-(|T^) and g n>q is the above exact 
solution of eq. (||) with G q = 2 M " ) . Formula (^) is valid 
for n = 1,2, ... and g = 2,3, .... 

Noteworthy, that small correlations between the local 
values of energy flux density e and amplification coeffi- 
cient M lead to a non-zero intcrmittency exponent /j. ni i 
(which otherwise is zero and which usually is assumed to 
be zero). It follows from the definition of M rewritten in 
the form e(R, y) = e(r, x)M(R, x, y). Space averaging of 
this relation gives 



e(l-M) = (e-e)(M-M) 



(20) 



where e does not depend on x (or y) due to the total flux 
conservation. If correlations between e and M were ab- 
sent, the right hand side of (|2Q ) would be zero. It would 
imply M = 1 and // n) i = 0, in agreement with eq.(|3j) 
that neglects considered correlations. Positive correla- 
tions between e and M, i.e., a positive right-hand side of 
(|20|), implies that M < 1, which means in turn ^ nj i < 0. 
According to g , the g-dependence of the intermittency 
exponents is given by 



Mn, 9 = In q [arsh {q/a n ) - arsh (l/a„)] + /i„,i <? 



(21) 



g < q n 



u 



3/7n 



l/2+(l/4+(a n7n /3) 2 ) 



1/2 



1/2 



(which corresponds to the dimensionality D n q = 3 + 

fi>n, g -qn'n,q = 3 -Jnq 2 /\/q 2 + <%, so that D n , qnM = 0). 
It can be verified now that these formulas (have not been 
used above except the case of infinite n, eq.([13|)) are well- 
compatible with the current theory. This allows one to 
extend the predicted [i n ,q for all positive q. For n = 1, 
an extra constrain of the kind (\L&) should hold: 



F{a lmM , ai, ) « («„,, = ^ nq ) . 

It allows one to evaluate the single unknown parameter 
/i2 in formula (19h for fj, n ,q- Specifically, minimization of 
the norm of relative deviation between right-hand sides 
of ( |19| ) and ( |2l| ) for n = 1, subject to the constrain (22), 
gives [ii rs 0.215 for the intermittent turbulence. 



• COMPARISON WITH EXPERIMENT: 

The deviation of the above theoretical fi n q from the 
experimental data of [B does not exceed 8 [i q , where 
5 w /X2/2 ~ 0.023 is the small parameter. The devia- 
tion contains regular oscillations in n, which sweep is of 
the order of S /i q (see Fig.l). This indicates that there is 
an additional small effect that has not been taken into 
account above. It consists in some oscillations in the 
probability of blow-up branching (i.e., hot spot splitting) 
as the function of n. Stronger branching makes larger 
the local dimensionality D nq and reduces the local am- 
plification M. The regions of stronger branching repeat 
with some period in n, which implies a correlation be- 
tween consequent M. Calculation of such a correlation 



should improve the accuracy of theory (and probably put 
it beyond the accuracy of currently available experimen- 
tal data). 

• SUMMARY: This work has far reaching theo- 
retical implications. In particular, it indicates that the 
intermittency problem can be solved by appropriate per- 
turbation methods. The reason is that the major corre- 
lations causing intermittency can be described explicitly, 
while the remaining correlations are small and can be 
treated as perturbations. More specifically, it is shown 
that exactly solvable model (||),(|5|) can be appropriately 
"dressed" (see eq.(|l9|) with the environment) to describe 
very well the experimentally measured intermittency ex- 
ponents fi Hiq at all n and q for which these values are cur- 
rently available. Noteworthy, that the equations derived 
for scale-dependent intermittency exponents [i. n ,q do not 
contain any adjusting parameter (that should be deter- 
mined from comparison with experiment), but do have 
both the Kolmogorov and the intermittent solutions. 
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